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DISK- JET CONNECTION IN THE RADIO GALAXY 3C 120 
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ABSTRACT 

We present the results of extensive multi- frequency monitoring of the radio galaxy 3C 120 between 
2002 and 2007 at X-ray (2-10 keV), optical (R and V band), and radio (14.5 and 37 GHz) wave 
bands, as well as imaging with the Very Long Baseline Array (VLBA) at 43 GHz. Over the 5 yr of 
observation, significant dips in the X-ray light curve are followed by ejections of bright superluminal 
knots in the VLBA images. Consistent with this, the X-ray flux and 37 GHz flux are anti-correlated 
with X-ray leading the radio variations. Furthermore, the total radiative output of a radio flare 
is related to the equivalent width of the corresponding X-ray dip. This implies that, in this radio 
galaxy, the radiative state of accretion disk plus corona system, where the X-rays are produced, has 
a direct effect on the events in the jet, where the radio emission originates. The X-ray power spectral 
density of 3C 120 shows a break, with steeper slope at shorter timescale and the break timescale is 
commensurate with the mass of the central black hole based on observations of Seyfert galaxies and 
black hole X-ray binaries. These findings provide support for the paradigm that black hole X-ray 
binaries and both radio-loud and radio-quiet active galactic nuclei are fundamentally similar systems, 
with characteristic time and size scales linearly proportional to the mass of the central black hole. 
The X-ray and optical variations are strongly correlated in 3C 120, which implies that the optical 
emission in this object arises from the same general region as the X-rays, i.e., in the accretion disk- 
corona system. We numerically model multi-wavelength light curves of 3C 120 from such a system 
with the optical-UV emission produced in the disk and the X-rays generated by scattering of thermal 
photons by hot electrons in the corona. From the comparison of the temporal properties of the model 
light curves to that of the observed variability, we constrain the physical size of the corona and the 
distances of the emitting regions from the central BH. In addition, we discus physical scenarios for 
the disk-jet connection that are consistent with our observations. 

Subject headings: galaxies: active — galaxies: individual (3C 120) — galaxies: jets — physical data 
and processes — X-rays: galaxies — radio continuum: galaxies 
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1. INTRODUCTION 

Stellar-mass black hole X-ray binaries (BHXRBs) and 
active galactic nuclei (AGNs) are both powered by accre- 
tion onto a black hole (BH). In many cases these systems 
emit radiation over se veral decades of frequency a nd pos- 
sess relativistic jets ([Mirabel fc Rodriguez! Il999) . The 
above similarities in the basic generation of energy and 
observational properties have led to the paradigm that 
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these two systems are fundamentally similar, with char- 
acteristic time and size scales linearly proportional to 
the mass of the central BH (-10 M for BHXRBs and 
10 6 to 10 9 M for AGNs). Although this paradigm has 
given rise to the expectation that we might test mod- 
els of AGNs with observations of BHXRBs, such an ap- 
proach is unjustified until detailed, possibly quantitative 
connections between BHXRB systems and AGNs become 
well-established. The comparison of BHXRBs and AGNs 
is complicated by the fact that a single AGN usually 
does not show the entire range of properties that we 
wish to compare. For example, Seyfert galaxies are the 
AGNs that most resemble BHXRBs, but their radio jets 
tend to be weak and non-relativistic (e.g. lUlvestad et all 
1999). On the other hand, in radio-loud AGNs with 
strong, highly variable nonthermal radiation (blazars), 
the Doppler beamed emission from the jet at most wave- 
lengths masks the thermal emission from the accretion 
disk and its nearby regions. 

One well-established property of BHXRBs is the 
connection between accretion state and events in 
the jet. In these objects, transitions to high- 
soft X-ray states are associated with the emergence 
of very bright featur es that proceed to propagate 
down the radio jet (iFender Belloni. fe Gallol 120041 ; 
iFender. Homan. fe Bellonil 12009ft . Similarly, an accre- 
tion disk-jet connect ion was established in AGNs by 
iMarscher et al.l (|2002f ). who reported a relationship be- 
tween X-ray and radio events in the radio galaxy 3C 120 
(redshift of 0.033). During — 3 yr of monitoring of this 
object, four dips in the X-ray flux, accompanied by spec- 
tral hardening, preceded the appearance of bright super- 
luminal knots in the radio jet. This Fanaroff-Riley class 
I radio galaxy has a prominent radio jet that displays 
stron g variability in flux and jet structure (|Gomez et al.l 
[200lHWalker et al.|[200ll ). The jet lies at an angle ~ 20° 
to our line of sight, significantly wider than is typically 
the case for blazars (|Jorstad et al.ll200"5l ) . At optical and 
X-ray frequencies, 3C 120 possesses properties similar 
to Seyfert galaxies and BHXRBs, e.g., a p rominent iron 
emission line at a rest energy o f 6.4 keV ([Grandi et al.l 
119971 iZdziarski fe Grandil l200l . Hence, most of the X- 
rays are produced in the immediate environs of the ac- 
cretion disk: the corona, a hot wind, or the base of the 
jet. Since the superluminal knots are disturbances prop- 
agating down the jet, a link between decreases in X-ray 
production and the emergence of new superluminal com- 
ponents demonstrates the existence of a strong disk-jet 
connection. 

In this paper we present the results of extensive multi- 
frequency monitoring of 3C 120 between 2002 and 2007 
at X-ray energies (2.4-10 keV), optical R and V bands, 
and radio frequencies 37 and 14.5 GHz, as well as imag- 
ing with the Very Long Baseline Array (VLBA) at 43 
GHz. We use these data to obtain quan titative support 
for th e disk-jet connection claimed by Marscher et alJ 
(2002). Times of "ejections" are defined as the extrap- 
olated time of coincidence of a moving knot with the 
position of the (presumed quasi - statio nary) core in the 
VLBA images. iSavolainen et al.l (|2002l ) have found that 
ejection times are generally associated with an increase 
in flux at 37 GHz. Therefore, the 37 GHz flux should be 
anti-correlated with the X-ray flux if decreases in X-ray 



production are linked to ejections. Here, we confirm the 
existence of this anti-correlation and test for a relation 
between the amplitude of the X-ray dips and that of the 
associated 37 GHz flares. 

Another similarity between BHXRBs and Seyfert 
galaxies lies in their X-ray power spectral densi- 
ties (PSDs). The PSD corresponds to the power in 
the variability of flux on different timescales. The 
X-ray PSDs of Seyfert galaxies c an be fit by piece- 
wise power laws with a bre a k (|Nowak et al.l 119991: 
Uttlev. McHardv. fe Papadakisl 120021: iMcHardv et af] 
20041: iMarkowitz et alJ 120031: iPounds et all 120011 : 
Edelson fe Nandral 119991 ). This is generally similar to 
that of BHXRBs although in many cases the latter 
have more than one break and sharp peaks due to 
quasi-periodic oscillations. We use our data set to verify 
that the X -ray PSD of 3C 120 also has a break, as 
inferred by iMarshall et al.l (|2009l ). who calculated the 
PSD from somewhat less extensive observations that 
partially overlap with ours. 

The bulk of the optical continuum from non-blazar 
AGNs is thou ght to emanate from the accretion disk 
(|Malkanl [T983). In order to test this, we cross-correlate 
the optical flux variations of 3C 120 with those at X- 
ray energies. However, the coefficients of raw correla- 
tions may be artificially low due to the uneven sam- 
pling present in the data. In light of this, we use simu- 
lated light curves, based on the underlying PSD, to esti- 
mate the significance of the derived coefficients of both 
the X-ray/optical correlation and X-ray/37 GHz anti- 
correlation found in the 3C 120 data. 

In §2 we present the observations and data reduction 
procedures, while in §3 we describe the power spectral 
analysis and its results. In §4 we cross-correlate the X- 
ray and 37 GHz light curves and quantitatively inves- 
tigate the relation between the X-ray dips, radio flares 
and ejections, as well as discuss the X-ray/optical corre- 
lation. In §5 we discuss and interpret the results, focusing 
on the implications regarding the accretion disk-jet con- 
nection and the source of the optical emission in 3C 120. 
§6 presents conclusions, including comparison of 3C 120 
with BHXRBs 

2. OBSERVATIONS AND DATA ANALYSIS 

Table Q] summarizes the intervals of monitoring at dif- 
ferent frequencies for each of the three wave bands in 
our data set. We term the entire light curve "monitor 
data" ; shorter segments of more intense monitoring are 
described below. 

The X-ray light curves are based on observations of 
3C 120 with the Rossi X-ray Timing Explorer (RXTE) 
from 2002 March to 2007 May. We observed 3C 120 with 
the RXTE PCA instrument with typical exposure times 
of 1-2 ks. For each exposure, we used routines from the 
X-ray data analysis software FTOOLS and the program 
XSPEC to calculate and subtract an X-ray background 
model from the data and to fit the X-ray spectrum from 
2.4 to 10 keV as a power law with low-energy photoelec- 
tric absorption by intervening gas in our Galaxy. For the 
latter, we use d a hydrogen column density of 1.23 x 10 21 
atoms cm" - 2 (lElvis. Lockman. fe Wilkeifl989h . 

The sampling of the X-ray flux varied. Normally, ob- 
servations were made 2-3 times per week except during 
8-week intervals each year when the radio galaxy is too 
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close to the Sun's celestial position to observe safely. In 
order to sample shorter-term variations, between 2006 
November and 2007 January we obtained, on average, 
four measurements per day for almost two months. We 
refer to these observations as the "medium" data. XMM- 
Newton observed 3C 120 quasi-continuously for about 
130 ks on 2002 December 13 and 14, during which all 
instruments were operating normally. The data were 
processed with the latest software (SAS version 5.3.3). 
Light curves were extracted in two energy bands, 0.3-10 
keV and 4-10 keV, and were background-subtracted and 
binned to 100 s time intervals. The 4-10 keV light curve, 
re-binned to an interval of 0.01 day, were used as the 
"longlook" data. We use the 4-10 keV data since the en- 
ergy range is similar to that of our RXTE data (2.4-10 
keV). Figure Q] presents these three data sets. 

We also monitored 3C 120 in the optical R and V 
bands over a portion of the time span of the X-ray 
observations. The majority of the measurements in R 
band are from the 2 m Liverpool Telescope (LT) at La 
Palma, Canary Islands, Spain, supplemented by obser- 
vations at the 1.8 m Perkins Telescope of Lowell Ob- 
servatory, Flagstaff, Arizona. The V-band photometry 
was obtained with the 0.4-m telescope of the Univer- 
sity of Nebraska. On each night a large number of one- 
minute images (~ 20) were taken and measured sep- 
arately. Details of the o bserving and reduction pro- 
cedure are as described in IKlime k. Gask ell. fc Hen drick 
(2004) . Comparison star mag nitudes were calibrated us- 
ing Doroshcnko ei~aT1 (l2?ffll . To minimize the effects 
of variations in the image quality, fluxes were measured 
through an aperture of 8 arcseconds radius. The errors 
given for each night are the uncertainties in the means. 
The Miller Observatory observations were taken with 
a 0.4-m telescope in Nebraska and reduced in similar 
manner as were the University of Nebraska observations. 
Observations at the Shanghai Astronomical Observatory 
were obtained with the 1.56-m telescope at Sheshan Sta- 
tion with standard Johnson-Cousins V, R, and I filters, 
and all the magnitudes were scaled to the V passband. 
The reductions were as for the Nebraska observations. 
Early observations (2002 to 2004) carried out at the 
Perkins Telescope were in V-band. For these m e asure - 
ments, we used stars D, E, and G from lAngiond (|1971l ) 
to calculate the V-magnitudes. We use three comparison 
stars in the field of 3C 120 to calculate the R-magnitudc. 
We determined the R-magnitudes of these three stars 
based on ~20 frames obtain ed within 2 yr. W e use the 
flux-magnitude calibration of lMead et al.l (|1990f ) and cor- 
rect for Galactic extinction for both R and V bands. 

Since the optical sampling in R or V band individually 
is not as frequent as the X-ray sampling, we construct 
a better-sampled optical light curve by combining these 
two bands. We find that the R and V band light curves 
have 38 data points that were measured within 0.5 day 
of each other. These data are shown in Figure O The 
equation of the best-fit line is Fr — 0.96i*V + 1-96, where 
Fy and Fr are the fluxes in V and R band, respectively. 
We convert the V band fluxes into R band using this 
equation. We also add some data from Lowell Observa- 
tory while constructing the combined light curve, which 
we present in Figure [3J 

We have compiled a 37 GHz light curve with data from 
the 13.7 m telescope at Metsahovi Radio Observatory 



Finland. The flux density scale is set by observations of 
DR 21. Sources 3C 84 and 3C 274 are used as secondary 
calibrators. A detailed description on the data reduction 
and analysis is given in Terasranta ct al. (1998). We also 
monitored 3C 120 at 14.5 GHz with the 26 m antenna of 
the University of Michigan Radio Astronomy Observa- 
tory. Details of the calibration an d analysis techniques 
are described in lAller et alj (|1985l ). At both frequencies 
the flux scale was set by observations of Cassiopeia A 
(jBaars et al.lll977D . 

Starting in 2001 May, we observed 3C 120 with the 
Very Long Baseline Array (VLBA) at 43 GHz at roughly 
monthly intervals, with some gaps of 2-4 months. The 
sequence of images from these data provide a dynamic 
view of the jet at an angular resolution ~0.1 milliarcsec- 
onds (mas) in the direction of the jet, corresponding to 
0.064 pc for Hubble constant H Q = 70 km s" 1 Mpc -1 . 
W e processed the data in the same manner as described 
in lJorstad et al.l (|2005D . We modeled the brightness dis- 
tribution at each epoch with multiple circular Gaussian 
components using t he task MODELFIT of the software 
package DIFMAP (|Shepherd et al.lll997l) . This repre- 
sents the jet emission downstream of the core by a se- 
quence of knots (also referred to as "components" ) , each 
characterized by its flux density, FWHM diameter, and 
position relative to the core. The apparent speeds of the 
moving components with well-determined motions are all 
within, 4.0c ± 0.2c. The ejection time To is the extrap- 
olated time of coincidence of a moving knot with the 
position of the (presumed stationary) core in the VLBA 
images. Table [5] lists the ejection times and Figure 2] and 
[5] display the VLBA images. Fig. [5] presents the X-ray, 
optical and radio light curves. In the top panel of the 
figure, the arrows represent the times of supcrluminal 
ejections, while the line segments perpendicular to the 
arrows show the uncertainties in the values of Tq. Fig- 
ure [7] displays the X-ray and optical light curves in the 
interval between 2006 October and 2007 April. It shows 
the best-sampled pair of dips in the X-ray and optical 
fluxes. 

3. POWER SPECTRAL ANALYSIS AND RESULTS 

We use a varia nt of the Power Spectrum Response 
method (PSRESP lUttlev. McHardv. fc Papadakiil2002T) 
to determine the intrinsic PSD of the X-ray light 
curve. Our re a lizatio n of PSRESP is described in 
IChatteriee et al.l (|2008D . PSRESP gives both the best- 
fit PSD model and a "success fraction" F succ (fraction 
of simulated light curves that successfully represent the 
observed light curve) that indicates the goodness of fit of 
the model. 

At first we fit a simple power-law model to the X-ray 
PSD, but found that the value of F succ was unacceptably 
low. This implies that a simple power-law is not the best 
model for this PSD. Then we fit a bending power-law 
model (broken power law with a smooth break) to the 
X-ray PSD, 

P(u) = Au aL [l + (— )(<« -«*)]-!. (i) 

Here, A is a normalization constant, vb is the break fre- 
quency, and an and ai are the slopes of the power- 
laws above and below the break frequency, respectively 
(jMcHardv et al.ll2004f h During the fitting, we varied vb 
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from 10~ 9 to 10~ 5 Hz in steps of 10 05 , «h from —1.5 
to —3.0 in steps of 0.1, and «l from —1.0 to —1.5 in 
steps of 0.1. These ranges include the values of a found 
in the light curves of BHXRBs in the high-soft state, 
for which cei, » —1 and an b etween —2 and —3 (e.g. 
iRemillard. fc McClintockl 120061 ) . This procedure yields 
a much higher success fraction than the simple power- 
law model. Based on the model with the highest success 
fraction, we obtain a best fit with the parameters ql = 

-1.3±°;?, a H = -2.518;!, and logioM = -5.0518:2 Hz. 
The success fraction for this fit is high, 0.9. The intrin- 
sic scatter of power is very high for slopes below —2 due 
to red noise leak. As a result, trial PSD models with 
a high frequency slope of any value between —2 and —3 
give very similar success fractions. This causes the larger 
uncertainty in the value of the high frequency slope rel- 
ative to that at lower frequencies. Figure [8] presents this 
best-fit model and the corresponding PSD. As seen in 
the figure, the high frequency part of the PSD is dom- 
inated by Poisson noise. That is because this part of 
the PSD is generated from the longlook light curve, and 
fluxes in the longlook light curve have larger uncertain- 
ties owing to shorter exposure times than those in the 
other light curves. The figure shows that when the esti- 
mated Poisson noise is added to the best-fit model PSD, 
it matches the observed PSD quite well. The nature of 
the X-ray PSDs and the values of break frequencies found 
previously in Seyfert galaxies ar e consistent with thos e 
of BHXRBs in high-soft state (jMcHardv et alJ 120051 ). 
This is consistent with accretion rates in all these ob- 
jects which are above 10%, more than the typical value 
of 2% indicating the transition from l ow to high state 
in galactic BHXRBs (jMaccaronel 120031 ). In 3C 120, the 
accretion rate ~30% which implies that it is also a high 
state system. 

The PSD break frequency in BHXRBs and Seyfert 
galaxies scales inversely with the mass of the black 
hole while being proportional to th e accretion rate 
(lUttley. McHardv. fc Papadakis) 120021 : iMcHardv et al l 
20041 120061 : lEdelson fc Nandral 119991 : iMarkowitz et all 
20031 ). Using the best-fit values and uncertainties in the 
relation between bre ak timescale, BH mas s, and accre- 
tion rate obtained by IMcHardv et all (|2006[ ), we estimate 
the expected value of the break frequency in the X-ray 
PSD of 3C 120 to be l(r 5 - 0±0 - 7 Hz. Hence, our derived 
break frequency lies within the expected range. Here 
we adopt a BH mass of 5.5 x 10 7 from e mission-line 
reverberation mapping ([Peterson et al.ll2004T ) and a bolo- 
metr ic luminosity of the big blue bump of 2.2 x 10 45 ergs 
s^ 1 (IWoo fc TJrrvl[2002T ). 

We note that the high and low frequ ency slopes, as 
well a s the break frequency calculated by Marshal l et al.l 
(2009) (-2.1 ±0.4, -1.2, io- 5 - 75±0 - 43 Hz), are consistent 
with the best-fit values in this paper within the uncer- 
tainties. The difference in the central values probably 
results from our use of a different longlook light curve, 
which is from the EPIC instrument onboard XMM- 
Newton. The sensitivity of this detector is higher than 
the PCA of RXTE, allowing our light curve to have a 
much better sampling rate. This is useful for extending 
the lower end of the frequency range of the PSD, leading 
to a more accurate measurement of the high frequency 
slope and break frequency. 



4. CROSS-CORRELATION ANALYSIS AND RESULTS 

We e mploy the discrete cro ss-correlation function 
(DCCF lEdelson fc Krolikl [l988h method to search for 
connections between variations at pairs of wave bands. 
We determine the significance of the correlations by 
performing the analysis on 100 simulated light curves 
generated randomly with a PSD matching that of the 
real data. Deta i ls of this procedure are described in 
iChatterjee et al.l ((2008) . In the present case, we corre- 
late simulated X-ray light curves with the observed op- 
tical/radio light curves to calculate the significance. 

4.1. X-ray/Radio Correlation 

The 37 GHz light curve has an average sampling fre- 
quency of about once per week. We bin the X-ray and the 
37 GHz light curves in 7-day intervals before performing 
the cross-correlation so that the light curves being com- 
pared have similar sampling frequency. As determined 
by the DCCF (top panel of Figure [H|), the X-ray flux 
variations are anti-correlated with those at 37 GHz in 
3C 120. The highest amplitude of the X-ray versus 37 
GHz DCCF is at a value of -0.68 ± 0.11, which corre- 
sponds to a 90% significance level. The time lag of the 
peak indicates that the X-ray lead the radio variations 
by 120 ± 30 days. We have used the FR-RSS technique 
proposed bv lPeterson et al.l (|1998f ) to calculate the mean 
value and uncertainty of the cross-correlation time lag. 
This procedure gives quantitative support to the trend 
that is apparent by inspection of the light curve, i.e., X- 
ray dips are followed by appearance of new superluminal 
knots and hence enhancement in the 37 GHz flux. 

The bottom panel of Figure [5] shows the X-ray/37 GHz 
DCCF but without the data after 2006 April, in order to 
exclude the two deepest X-ray dips and highest ampli- 
tude flare at 37 GHz toward the end of our monitoring 
program. The X-ray/37 GHz anti-correlation remains, 
although the minimum of the DCCF is smaller, —0.4, 
with a significance level of 72%. The X-ray variations 
lead the radio by 80 ±30 days. Thus, the X-ray/37 GHz 
anti-correlation is not just the consequence of a singular 
event. The longer time delay of the major event in late 
2006 was caused by the longer time between the start of 
the outburst and the peak relative to other radio flares, 
especially the one at 2003.7. This is naturally explained 
as a consequence of higher optical depth of stronger radio 
outbursts. 

4.2. X-ray Dips and Radio Flares 

In order to determine the physical link between the 
accretion disk and jet, we check whether the amplitudes 
of the X-ray dips and associated radio flares are related. 
To test this, we calculate the equivalent width of each 
X-ray dip and the area under the curve of each 37 GHz 
flare to measure the total energy involved in the events. 
We approximate that the radio light curve is a super- 
position of a constant baseline of 1.5 Jy from the more 
extended jet and long-term flares. The baseline is cho- 
sen as approximately t he lowest level of the 3 7 GHz light 
curve. We then follow IValtaoj a et alJ (|1999l ) by decom- 
posing the baseline-subtracted light curve into individual 
flares, each with exponential rise and decay. We use four 
parameters to describe each flare: the rise and decay 
timescales, and the height and epoch of the peak. Be- 
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fore the decomposition, we smooth the light curve using 
a Gaussian function with a 10-day FWHM smoothing 
time. The details of the decompo sition procedure are 
described in lChatterjee et all (|2008l ). 

The X-ray light curve has a long-term trend, i.e., the 
baseline is not constant. We define the baseline X-ray 
flux as a cubic-spline fit of the annual mean plus one 
standard deviation. Although this is an arbitrary defini- 
tion, this baseline reproduces reasonably well the mean 
flux level in between obvious dips and flares. There is a 
long-term trend in the spectral index variations as well. 
The X-ray spectral index a x , defined by f x oc v a * , where 
f x is the X-ray flux density and v is the frequency, varied 
between —0.5 and —1.1, with an average value of —0.83 
and a standard deviation of 0.10 over the 5.2 yr of obser- 
vation. In addition to the short term fluctuations, there 
is a long-term trend of increasing values of a x during this 
interval. We calculate a baseline in the same manner as 
for the X-ray flux variations to highlight the change in 
spectral index during the dips (see Figure fT0|> . Before we 
perform the above analysis to calculate the baseline, we 
smooth the X-ray light curve using a Gaussian function 
with a 10-day FWHM smoothing time. 

We proceed by subtracting the baseline and then mod- 
eling the dips as inverted exponential flares. Figure [TT] 
shows the model fits to both the X-ray and 37 GHz light 
curves. Table[2]lists the parameters of the X-ray dips and 
the corresponding 37 GHz flares along with the times of 
superluminal ejections. Since the correlation functions 
with and without the data in 2007 (during the large am- 
plitude 37 GHz flare) show that the X-ray dips lead the 
37 GHz flares by 120 ± 30 and 80 ± 30 days, respec- 
tively, and we have identified 15 significant dips during 
the 5.2 yr of monitoring, we assume that a radio flare 
that peaks between 20 and 180 days after an X-ray min- 
imum is physically related to it. Using the results of 
this analysis, we plot the equivalent width of X-ray dips 
versus the energy output of the corresponding 37 GHz 
flares in Figure [TJ The 37 GHz flare at 2003.72 has 
a wide decaying wing and there are two X-ray dips (at 
2003.58 and 2003.82) that may be related to this flare 
as well as another smaller 37 GHz flare (at 2003.92) on 
top of that wing. Similarly, there are two X-ray dips at 
2005.12 and 2005.39, and there are four corresponding 37 
GHz flares very close to each other in time at 2005.23, 
2005.36, 2005.57, 2005.80. In the above cases, we plot- 
ted the total energy output of the group of 37 GHz flares 
against the total equivalent width of the corresponding 
group of X-ray dips. The energy output of the flares and 
the dips were corrected by adding the residuals shown in 
Figure [TT1 The plot in Figure [T2l reveals a positive corre- 
lation between the flare and dip strengths. The level of 
change in the accretion disk and/or corona is therefore 
closely related to the amount of excess energy injected 
into the jet. 

It can be seen from Table [5] that 14 out of 15 sig- 
nificant X-ray dips are followed by a superluminal ejec- 
tion. The time delay between the start of the X-ray dips 
and the time of "ejection" of the corresponding super- 
luminal knot varies between 0.03 yr to 0.46 yr, with a 
mean value of 0.18±0.14 yr. We plot the times of ejec- 
tion of new knots along with the corresponding times of 
start of the X-ray dips in Figure [13] The uncertainties 
in the start time are proportional to the widths of the 



dips and those in the ejection times are as given in Ta- 
ble [2] A straight line fits the data extremely well with 
small scatter, which indicates that there is a clear asso- 
ciation between X-ray dips and superluminal ejections. 
The best-fit line through the points has a slope of 1 and 
y- intercept of 0.19, which is consistent with the mean 
delay of0.18±0.14 y r given above. This s trongly sup- 
ports the proposition (Ma rscher et al.|[2002f ) that in 3C 
120, a decrease in the X-ray production is linked with 
increased speed in the jet flow, causing a shock front to 
subsequently move downstream. Generally, there is a 
close correspondence of superluminal ejections with 37 
GHz flares, but sometimes the decrease in flux farther 
out in the jet offsets much of the flux increase from the 
appearance of a new knot. In such cases, the correspond- 
ing increase of flux was not large enough to be detected 
in our decomposition of the smoothed light curve. This 
causes the minimum in the X-ray DCCF to be less sig- 
nificant than it would if old knots were to completely 
fade before new ones appear. Despite this complication, 
the X-ray/37 GHz cross-correlation (§4.1) verifies with 
an objective statistical method that radio events in the 
jet are indeed associated with X-ray dips. The apparent 
speeds of the moving components with well-determined 
motions are all similar, 4.0c ± 0.2c. Therefore, a knot 
moves a distance of 0.22 pc in 0.18 yr, projected on the 
plane of the sky. Since the angle of the jet axis of 3C 120 
to the line of sight ~20°, the actual distance traveled by 
the knot ~ 0.5 pc. Hence, we derive a distance ~ 0.5 
pc from the corona (where the X-rays are produced) to 
the VLBA 43 GHz core region. This confirms that the 
core is offset from the position of the BH. This is one of 
the few cases where we are able to probe the upstream 
region of the core using the time variable emission at a 
combination of wave bands. 

4.3. X-ray/Optical Correlation 

We bin the X-ray and the combined optical light curves 
in 2-day intervals before performing the cross-correlation, 
so that the light curves being compared have similar 
sampling frequency. As determined by the DCCF (Fig- 
ure I14|) . we find that the X-ray variations are very 
strongly correlated with those at optical wavelengths 
in 3C 120. The peak X-ray versus optical DCCF is 
0.80±0.07, which corresponds to a 99% significance level. 
The position of the peak of the correlation function in- 
dicates the relative time delay between the variations at 
the two wavelengths. In this case, the peak is very wide, 
so that the value of the relative time delay can not be 
easily estimated from the DCCF plot. We used t he FR- 
RSS technique proposed bv lPeterson et al.l (|1998l ) to cal- 
culate the mean value and the uncertainty of the cross- 
correlation time lag. This method indicates that the X- 
ray variations lead the optical by 0.5±4 days. The highly 
significant correlation and short time delay between the 
X-ray and optical variations indicates that emission at 
these wave bands is at least partially co-spatial. The 
asymmetry in the correlation function at time delays 
above ~100 days indicate that the X-ray variations lead 
those in the optical at these longer timescales. 

To characterize the variation of the X-ray/optical time 
lag over the years, we divide both light curves into two 
intervals, 2004 July to 2005 May and 2005 June to 2007 
May, and repeat the DCCF analysis on each segment. 
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The result (Figure fTS]) indicates that in the first segment, 
the X-ray variations lead those in the optical by ~25 days 
while in the second segment the correlation function is 
similar to what we obtained for the entire time interval, 
with similar non-zero time delay (Figure [T4")) . This vari- 
ation of the time lag over the years may be the cause of 
the observed wide peak in the correlation function. 

Between 2006 November and 2007 January, the X-ray 
light curve ("medium" light curve described in §2) was 
sampled 4 times per day and the combined optical light 
curve has a sampling rate of twice per day, on average. 
We bin these light curves in 0.5 day intervals. We cross- 
correlate the binned light curves in order to compare 
the correlation function with that of the long-term light 
curves. The correlation function, shown in Figure fTBl has 
a similarly significant correlation coefficient and time de- 
lay. The similar values of correlation coefficient and time 
delay using very well-sampled light curves illustrates the 
robustness of the correlation result. 

The X-ray to op tical time delay calculated by 
iMarshall et al.l (|2009f ) has the same sign as determined 
in this paper, although the magnitude is significantly dif- 
ferent (28.73j£g? days). We note that (1) This work 
includes more data points in the R band light curve, in- 
cluding the period 2007.0-2007.36, and 2) We construct 
a combined light curve using the R band data as well as 
scaled V band data. We use this more completely sam- 
pled light curve for the correlation analysis. This, in ad- 
dition to the longer time lag from 2004.5 to 2005.5, can 
account for the difference betwee n our result and that 
obtained by IMarshall et all (|2009f ). 

The interpolated cross-correlation function (ICCF) is 
another method to calculate the cross-correlation of un - 
evenly sampled discrete data (jGaskell fc Petersonlll987h . 
We have repeated all the above correlation analysis us- 
ing the ICCF and have found that the DCCF and ICCF 
results are consistent with each other for the data pre- 
sented in this work. 

We have cross-correlated the hard (4-10 keV) and soft 
(0.2-4 keV) longlook X-ray light curves, finding an ex- 
cellent correlation with nearly zero time lag (Figure I17j) . 
This indicates that in 3C 120, variations in soft and hard 
X-rays follow each other very closely. 

5. DISCUSSION 

5.1. Disk- Jet Connection 

The physical cause of the connection between events 
in the central engine and the jet of BHXRBs and AGNs 
is currently a matter of considerable speculation. If 
the jet is magnetically lau nched from the accretion disk 
(Blandfor d fc Pavnd 119821 ) . then there must be a link 
between the magnetic state at the base of the jet and 
the accretion state in the inner disk. One scenario, pro- 
posed for BHXRB s bv lLivio. Pringle. fc Kingl (|2003l) and 
iKing et ahl (|2004T ). involves a change in the magnetic 
field configuration in the inner disk from a turbulent 
condition in the high-soft state (when the X-ray flux is 
relatively higher and softer) to mainly poloidal in the 
low-hard X-ray state. The turbulence is needed for vis- 
cous heating, which in the BHXRB case leads directly to 
bright X-ray emission with a soft spectrum. In an AGN 
it causes strong ultraviolet emission, which is Compton 
scattered in the corona to a hard X-ray spectrum. If 



the field switches from chaoti c to mainly poloidal, which 
iLivio. Pringle. fc Kind (|2003[ ) suggest can occur by ran- 
dom episodes of near-alignment of the field in the rela- 
tively small number of turbulent cells in the inner disk, 
then the radiation in the inner disk will be quenched 
at the same time as energy flow into the jet is promoted. 
The transition back to the turbulent, radiative inner disk 
of the high-soft state would need to involve a surge of en- 
ergy injected into the jet in order to send a shock wave 
(|Miller- Jones et al~ll2005f ) down the jet. Perhaps global 
magnetic reconnection could cause this, but no detailed 
MHD model has been published to date. 

Alternatively, it is possible that the "corona," where 
the X-ray emission seen in AGN supposedly arises from 
Compton up-scattering of so fter accretion-disk photons, 
migh t be the base of the jet (|Markoff. Nowak. fc Wilms 
2005). If this is the case, then the X-ray flux will be re- 
lated to the number of electrons residing there and avail- 
able for scattering to create X-rays. The mass loading of 
the jet should also affect the asymptotic Lorentz factor of 
the flow downstream if th e jet is magnetically driven (e.g. 
IVlahakis fc K onigl 200l. The same decrease in electron 
number that causes a drop in scattered X-ray emission 
near the disk would lead to a time-delayed increase in the 
speed of the jet downstream. The flatter-spectrum non- 
thermal X-ray emission from the downstream jet would 
then play a larger relative role in the X-ray emission, 
causing the observed hardening of the spectrum during 
the dips. The increase in flow speed of the jet could form 
a shock wave, seen as a superluminal radio knot. It is dif- 
ficult to speculate why the mass loading should change, 
since we do not understand the processes by which ma- 
terial from the disk and/or ergosphere are injected into 
the base of the jet. However, observations of the mi- 
croquasar GRS 1915+105 suggest that outflow of matter 
from the disk switch es from mainly a wind to the jets 
(jNeilsen fc Ledl2009l ). In an AGN with a magnetically 
driven jet, it may be the case that lower mass injection 
into the jet actually enhances the jet emission owing to 
an increase in the flow speed driving a shock wave down 
the jet, as described above. 

5.2. Source of Optical Emission 

The strong correlation between the X-ray and opti- 
cal variations in 3C 120 implies that the emission at 
both wave bands arises from the same general region. 
Since the X-rays originate in the corona, the optical 
emission i s probably th ermal emission from the accre- 
tion disk (Malkan 1983). We can reject the alternative 
hypothesis that the main component of optical emission 
is synchrotron radiation from the jet. In that case, the 
emission should be significantly polarized, contrary to 
observation s showing the optic al linear polarization to 
be < 0.3% (porstad et al.l 120071) . Furthermore, the anti- 
correlation of the X-ray and 37 GHz emission (the latter 
of which is produced in the jet) implies that any opti- 
cal synchrotron emission from the jet should also anti- 
correlate with the X-ray flux, contrary to the strong ob- 
served correlation. 

Although the X-rays are predominantly produced by 
inverse Compton (IC) scattering of the thermal opti- 
cal/UV seed photons from the accretion disk by hot, 
but non-relativistic, electrons in the corona, the opti- 
cal/UV emission and the X-rays are tied together by 



7 



another mechanism: some of the optical/UV radiation 
is produced by heating of the accretion disk by X-rays 
produced in the above process ("feedback" mechanism). 
The amplitude of the short timescale (~days) variations 
in the X-ray light curve is larger than that in the optical. 
This can happen if a significant fraction of the optical 
emission is due to heating of the disk by X-rays, since 
this reprocessing may smooth out the short timescale 
variations. Feedback can make the time delay more diffi- 
cult to define by producing a fraction of the optical/UV 
photons with different temporal properties from that of 
the direct emission. 

There is a weak correlation between the X-ray spectral 
index and the X-ray flux, indicating that the spectrum 
becomes harder during a decrease in the production of X- 
rays fFiglTOl). T his ki nd of "pivoting" has been seen be- 
fore (jOele et al.ll2005t iMaraschi et alJll991f >. The X-rays 
may be produced mainly by up-scattering of UV photons 
and not optical photons, co nsistent with the correlation 
found by lOgle et alJ (|2005l ) between the X-ray and UV 
variations in 3C 120. This could occur if the flux of op- 
tical photons reaching the corona is much smaller than 
that of UV photons. Such a scenario is likely if the corona 
is small such that the region where the UV photons are 
generated is much closer to the corona than the region 
where the majority of the optical photons are produced. 
Any disturbance propagating outwards in the accretion 
disk will cause a change in the UV flux (and a resultant 
nearly immediate change in the X-ray emission) followed 
by a similar change in the optical flux. The sign of the 
time delay will switch if the disturbance propagates in- 
wards. The observed mean time delay between X-ray and 
optical variations may be due to such propagation time 
delays. If we adopt typical para meters such as ~ .1 for 
accretion efficiency, L/L e d ~ 0.3lOgle et al.l (|2005| ). and 
L bol ~ 2x 10 45 ergs/sec (|Woo fc Urry||2002[ ) and neglect 
General Relativistic effects, the region in the accretion 
disk where the emission peaks in the extreme UV/soft 
X-ray range (A = 10 nm) should be very close to the in- 
nermost stable orbit (~5r 9 , where r g is the gravitational 
radius of the black hole). The region where the emission 
peaks in the optical (A = 600 nm) lies at ~ 1000r g . For 
3C 120, with a BH mass of - 5xl0 7 M Q , 1000 r g is 
equivalent to ~4 light days. We consider a model that 
includes both time delays and coupling of emission from 
the corona and disk. 

We have performed a theoretical calculation to pro- 
duce multi-wavelength light curves from a disk-corona 
system. We assume that (l) the temperature of the disk 
changes with its radial distance from the center accord- 
ing to T~ r^' A (jShakura fc Sunvaevl Il973h . (2) each 
annulus radiates as a perfect black-body, (3) the X-rays 
are produced in a spherical region (corona) close to the 
center of the disk by inverse-Compton scattering of the 
disk photons that reach the corona, and (4) lower energy 
radiation (UV-optical) is produced both by blackbody 
radiation in the disk and reheating of the disk by the 
X-rays from the corona. For simplicity, we assume that 
the energy distribution of electrons is uniform through- 
out the volume of the corona. The fraction of the UV- 
optical photons that reach the corona are up-scattered 
to X-rays and a part of the X-rays re-heat the disk to 
produce additional thermal photons. We treat the X-ray 



production in the corona as a reflection that increases the 
energy of the radiation. The computer code that we use 
for the model also follows the travel time of photons from 
different radii of the disk to the corona and vice versa. 
The temporal nature of the emission from the corona and 
the part of the disk close to it are very similar owing to 
the small light travel time. Figure [TS] shows the varia- 
tion of the total UV and optical intensity emitted by an 
annulus at radius r of thickness Ar oc r, i.e., B v r ver- 
sus r. This shows that the annulus that produces the 
largest amount of UV radiation ("region UV" ) lies ~5 r g 
from the center, and for optical emission ('region Op") 
it is at ~75 r g . This distance is much smaller than what 
would be expected if most of the optical radiation were to 
come from radii near where the Planck curve B v peaks 
at R band. These latter radii lie ~4 light-days (1000 
r g ) from the center, while regions close to the center are 
at higher temperature, and therefore emit at a higher 
blackbody intensity at all wavelengths. Although B v in- 
creases monotonically with decreasing radii, B v r peaks 
at a radius ~75 r g at R band. 

We use a computer code to model this disk-corona sys- 
tem and then introduce a disturbance in the temperature 
of the accretion disk that propagates at a speed ~0.1c 
from the center toward the outside or vice versa. Due 
to the disturbance, the temperature at a given annulus 
varies with time according to a Gaussian profile. This 
causes a flare in the emission of the entire system at all 
wavelengths, although the flare starts and peaks at dif- 
ferent times at different wavelengths. We produce these 
light curves at X-ray, UV, and optical bands, includ- 
ing delays from internal light-travel time. The nature 
of the flare and relation between flares at multiple wave 
bands depend on the speed and direction of propagation 
of the disturbance since emission from different annuli 
of the disk have different wavelength dependence. The 
top panel of Figure [TH shows the X-ray and optical light 
curves from the above calculation when the disturbance 
propagates toward smaller radii. As a result, the emis- 
sion from the outer disk flares before that closer to the 
BH. Therefore, the optical variation leads that in the 
X-rays. In this simulation, the time delay between the 
peaks of the X-ray and optical light curves ~4 days. 

In 3C 120, the time delay is centered on 0.5 day with an 
uncertainty of ±4 days, hence the direction of the time 
delay cannot be specified with certainty. The correlation 
function in Figure [16] also shows a relatively broad peak 
centered on zero, with a similar uncertainty in the time 
delay. The cross-correlation function for longterm light 
curves sometimes show a broad peak due to changes in 
the direction and magnitude of the time delay over the 
timespan of the light curves. If the broad peak in these 
correlation functions is due to the above-mentioned rea- 
son, it may indicate a dichotomy in the speed of propa- 
gation of the disturbances and their directions, i.e., the 
uncertainty corresponds to an actual range of time de- 
lays. This also constrains the size of the corona. For 
example, if the corona were spread such that the regions 
Op and UV were at the same distance from the corona, 
then these regions would contribute equally to the flux 
of seed photons that are up-scattered to X-rays. In that 
case, the X-ray flares would be much broader than ob- 
served and no optical/X-ray time delay would be present. 
This is contrary to the observation of relatively sharp X- 



ray flares. In addition, although overall we find a short 
mean time delay of 0.5 day (i.e., if the uncertainty does 
not correspond to an actual range of time delays but 
rather is consistent with zero), in some of the individ- 
ual flares the time delay is longer (~25 days; see top 
panel of Figure [15]). Also, if optical and UV photons 
were up-scattered to soft and hard X-rays, respectively, 
a significant "soft-hard" time delay would be expected. 
But from Figure [I?] it is clear that the soft-hard time 
lag cannot be more than 0.25 day, which is essentially 
unobservable with the longterm light curves. From the 
above discussion, we conclude that the corona is situated 
close to the center of the accretion disk, and its size is 
such that the flux of photons reaching the corona from 
region Op is negligible with respect to that from region 
UV. Based on the solid angle subtended, the coronal ra- 
dius should be less than 40r 9 for the flux from region Op 
to be less than 10% of that from region UV. 

In Figure [ini the two panels show results for different 
levels of feedback (the fraction of the X-rays produced 
in the corona that reflect toward and re-heat the disk). 
The top panel has feedback fraction equal to zero, i.e., 
none of the X-rays re-heat the disk. We can see from the 
figure that, as feedback is introduced, the resulting time 
delay becomes smaller and less precise. Figure I2TJ1 shows 
the light curves when the disturbance is propagating out- 
wards from the BH. In this expected, the X-ray 
variations lead those in the optical (by similar amount as 
in Figure [T9|) for the same propagation speed. The panels 
show, as above, that including feedback makes the time 
delay shorter. Hence, feedback from the corona may also 
contribute to the range of time delays that we see in the 
data. 

In another possible scenario, the X-ray variability is 
caused by intrinsic changes in the hot electrons in the 
corona and the UV-optical changes are due to feedback, 
i.e., there is no intrinsic variability in the accretion disk. 
In this case, the X-ray/optical time delay will solely be 
due to light travel time from the corona to the accretion 
disk. This will produce time lags o f a small fraction of 
a day (jKazanas fc Navakshinl 120011 ). which is too small 
to observe with the sampling of the data used in this 
work. The observed essentially zero time lag between the 
long-term X-ray and optical light curves (0.5 ± 4 days) 
and smaller variability amplitude in the short-term op- 
tical light curve than that in the X-rays are consistent 
with the above scenario. In fact the steepening of the X- 
ray spectrum and decrease in the mean X-ray flux level 
after 2006.0 (see Fig. [Til]) suggest a long-term steepen- 
ing of the energy distribution of electrons in the corona. 
But the longer X-ray/optical time delay over a signifi- 
cant portion of the dataset (Figure [15]) and comparable 
X-ray/optical B-b and variability amplitud e in the long- 
term light curves (jDoro shcnko ct al. 2009) with B-band 
variations sometimes leading, indicate that propagation 
of disturbances in the accretion disk must produce at 
least part of the X-ray/optical variability in this radio 
galaxy. 

Figure I2T1 shows a sketch of the accretion disk-corona 
system as derived in this paper. The black filled cir- 
cle is the position of the BH, the temperature of the 
accretion disk is shown by gray scale with lighter color 
meaning higher temperature, and the larger circular area 
filled with dots is the corona. "Region UV" and "Region 



Op" are shown as thick solid lines on the accretion disk. 
The radius of the corona and the distance of the relevant 
emission regions from the BH are shown in units of r g . 

From the above model and the correlation results, we 
conclude that most, if not all, of the optical emission in 
3C 120 is produced in the accretion disk. The X-rays 
are produced by scattering of (mostly) UV radiation in 
the corona. Optical/UV emission due to re- heating of 
the disk by the X-rays is also a possible ongoing mecha- 
nism that may cause the spread of the time delay. The 
presence of sharp X-ray flares and the X-ray/optical time 
delay of weeks in some individual flares indicate that the 
corona must lie sufficiently close to the BH that the flux 
of seed photons reaching the corona is dominated by UV 
light. 

The light curve during the dip at X-ray and optical 
wave bands between 2006 October and 2007 April, dis- 
played in FigurefT] is very well-sampled. It is evident that 
the optical flux starts to decrease after ~MJD 4025, ~ 40 
days earlier than the decrease starts in the X-rays. The 
minimum in optical occurs ~40 days earlier than that in 
the X-rays as well. This indicates that the decrease was 
caused by a disturbance propagating from the outer radii 
of the accretion disk toward the BH. If this disturbance is 
a thermal fluctuation propagating inward, then it should 
have an effective speed < 0.01c to cause a time delay of 
^40 days. This is one order of magnitude higher than the 
sound speed for a gas pressure dominated disk of temper- 
ature < 10 5 K. Therefore, the ab ove disturbance c annot 
be transmitted by sound waves (Kro lik et al.|[l99lh un- 
less the relevant regions of the disk are dominated by 
radiation pressure, which cannot be ruled out in 3C 120 
given its high accretion rate, nearly the Eddington value 
(0.3 Le)- On the other hand, t his scenario inconsis- 
tent with the model proposed by IKing et al.l (|2004f) . in 
which large scale alignment of poloidal magnetic field 
in the inner accretion disk from random fluctuations 
causes the decrease in the X-ray flux. Such alignment 
occurs at a timescale 2 R / H k{R 3 / GM) 1 / 2 , where R/H 
is the radi us to thickness ratio of the disk, k is a con- 
stant ~10 (|Tout fc Pringldll992fc IStone et al.lll996D . and 
(R 3 /GM) 1 / 2 is the disk dynamical timescale at radius R. 
This alignment timescale, for a 10M Q BH, is few seconds, 
which translates to ~50 d ays for a ~ 5 xAJFMq BH in 3C 
120. Model light curves of lKing et~aTl (|2004D also contain 
short timescale, small amplitude fluctuations on top of 
the big flares and dips caused by small-scale alignment 
of the poloidal magnetic field, similar to that observed in 
the light curves of 3C 120 presented in this paper. 

6. SUMMARY AND CONCLUSIONS 

This paper presents well-sampled, 5-yr-long light 
curves of 3C 120 between 2002 and 2007 at X-ray, op- 
tical, and radio wavebands, as well as monthly images 
obtained with the VLB A at 43 GHz. We have calculated 
the broad-band X-ray power spectral density (PSD) of 
3C 120. The X-ray and radio light curves were decom- 
posed into individual dips and flares, respectively, and 
the properties of the contemporaneous dips and flares 
were compared. That comparison, as well as, the X- 
ray/radio correlation analysis were used to find a con- 
nection between events in the jet and the accretion disk. 
We have also cross-correlated the X-ray and optical light 
curves and have determined the significance of the cor- 
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relations with simulated light curves based on the PSD. 
Using a numerical model, we have produced X-ray and 
optical light curves from the accretion disk-corona sys- 
tem. From the comparison of the temporal properties of 
the simulated light curves to that of the observed vari- 
ability, we infer the physical size of the corona and the 
distances of the emitting regions from the central BH. 
Our main conclusions are as follows: 

(1) The X-ray PSD of 3C 120 is best fit by a bending 
power-law model where there is a smooth change in the 
slope above a break frequency. The best fit value of break 
frequency for 3C 120 is 10~ 505 Hz, which agrees very well 
with the relation between break timescale, BH ma ss and 
accretion rate obtained by McHardv et alj (|2006f ) span- 
ning a range of BH mass from 10 M to 10 9 M & . This 
indicates that the accretion process in 3C 120 is similar 
to that of the BHXRBs. 

(2) The X-ray and 37 GHz variations are anti-correlated 
with the X-ray leading the radio by ~120 days. The 
anti-correlation remains even if the large X-ray dips and 
radio flares after 2006 are excluded. 

(3) Almost all X-ray dips are followed by the ejection 
of a new knot in the VLBA images. This and the anti- 
correlation mentioned in (2) imply that decrease in the 
X-ray production is linked with increased speed in the 
jet flow, causing a shock front to form and move down- 
stream. This property of 3C 120 is also similar to the 
galactic black hole systems where transitions to high-soft 
X-ray states are associated with the emergence of very 
bright features that proceed to propagate down the radio 
jet. 

(4) We derive a distance ~0.5 pc from the corona (where 
the X-rays are produced) to the VLBA 43 GHz core re- 
gion using the average time delay between the start of 
the X-ray dips and the time of "ejection" of the corre- 
sponding super luminal knots. 

(5) The X-ray and optical vari ations in 3C 1 20 ar e 
very strongly corr e lated (see also IMarshall et all {2009); 
IDoroshenko et akl (|2009D ). This correlation, absence of 
significant optical polarization, and anti-correlation of X- 



ray and 37 GHz variations (the latter of which are pro- 
duced in the jet), together imply that the optical emis- 
sion is blackbody radiation from the accretion disk. 
(6) The X-rays are produced by scattering of (mostly) 
UV radiation in the corona. Comparison of simulated 
light curves from a disk-corona system and the observed 
variation imply that the corona must lie sufficiently close 
to the BH that the flux of seed photons reaching the 
corona is dominated by UV light. 
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Fig. 1. — X-ray light curves with different sampling rates. 
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Fig. 2. — Filled circles show the 38 data points that are measured within 0.5 day of each other in V and R band, along with the respective 
uncertainties. The dashed line represents the best fit straight line through these points, which is used for the V band to R band flux 
conversion. 
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Fig. 3. — Light curve of 3C 120 constructed by combining the V and R band light curves. 
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Fig. 4. — VLBI images of 3C 120 at 7 mm obtained from 2002 to 2004. The contours and color show the total and polarized intensity, 
respectively. The images are convolved with an elliptical Gaussian beam of size 0.36x0.15 mas at PA = —6°. The global peak over all 
maps is 1.52 Jy/Beam. The contour levels are 0.25, 0.5, 1.0, 64.0, 90.0% of the global peak. Individual moving knots are marked at 
selected epochs. 
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15 




Fig. 6. — Variation of X-ray flux, optical flux density and radio flux density of 3C 120 from 2002 to 2007. In the top panel, the arrows 
show the times of superluminal ejections and the line segments perpendicular to the arrows indicate the uncertainties in the times. 
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Fig. 7. — Points represent the X-ray and optical light curves between MJD 3950 and 4200 (2006 October and 2007 April) when the time 
sampling was dense during a minimum in the light curves. The curves represent the same data smoothed with a Gaussian function with a 
3 day FWHM smoothing time. 
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Fig. 8. — Result of application of the PSRESP method to the light curves. The PSD of the observed data at high, medium and low 
frequencies is given by the solid, dashed and dotted jagged lines, respectively, while the underlying power-law model is given by the thicker 
solid bent line. Points with error bars (open squares, solid circles and asterisks for high, medium and low frequency range, respectively) 
correspond to the mean value of the PSD simulated from the underlying power-law model (see text). The errorbars are the standard 
deviation of the distribution of simulated PSDs. The broadband power spectral density is best described by a bending power law with low 
frequency slope —1.3, high frequency slope —2.5 and break frequency 10 — 5 05 Hz. 
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Fig. 9. — Discrete cross-correlation function (DCCF) of the X-ray and radio monitor data. The time delay is defined as positive if the 
X-ray variations lag those at radio frequency. Top panel shows the correlation function for the entire data set and the bottom panel shows 
the same excluding the data during the major flare at 37 GHz in 2006-07 and the corresponding deep dip at X-ray energies. 
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Fig. 10. — Variation of X-ray flux and the X-ray spectral index during 2002-2007. Curves are the spline interpolation of the yearly mean 
of the data plus one standard deviation. In case of X-rays, variations (both positive and negative) are relative to the baseline. 
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Fig. 11. — X-ray and 37 GHz light curves. Curves correspond to su mme d flux after modeling the light curve as a superposition of many 
individual X-ray dips or radio flares and a baseline as shown in Figure flOl 
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Fig. 12. — Total energy output of the 37 GHz flares versus equivalent width of the corresponding X-ray dips. 
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Fig. 13. — Times of ejection of new VLB A knots versus times of start of X-ray dips from Table [2] The dashed line is the best-fit straight 
line through the points. The y-intercept of this line indicates the value of mean time delay between the start of dips and times of ejections. 
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Fig. 14. — Discrete cross-correlation function (DCCF) of the optical and X-ray monitor data for the entire 5 yr interval. The time delay 
is defined as positive if the variations at the higher frequency waveband lag those at the lower frequency. 
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Fig. 15. — Variation of the X-ray/optical correlation function across two intervals. The X-ray variations lead those in the optical by 25 
days during the first interval and the two variations are almost simultaneous during the second interval. 
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Fig. 16. — Discrete cross-correlation function (DCCF) of the optical and X-ray data between 2006 November and 2007 January, binned 
to an interval of 0.5 day. 




Fig. 17. — Discrete cross-correlation function (DCCF) of the soft and hard longlook X-ray data. The time delay is defined as positive if 
the variations at the hard wave band lag those at the soft wave band. 
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Fig. 18. — Total intensity of radiation coming from different annuli of an accretion disk versus the annular radii. The solid and dashed 
curves correspond to UV and optical wavelengths, respectively. 
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Fig. 19. — Simulated X-ray and optical light curves from an accretion disk-corona system using our theoretical calculation. The disturbance 
is propagating inwards towards the BH. The top panel corresponds to zero feedback and the bottom panel shows the same after including 
feedback. 
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Fig. 20. — Simulated X-ray and optical light curves from an accretion disk-corona system using our theoretical calculation. The disturbance 
is propagating outwards away from the BH. The top panel corresponds to zero feedback and the bottom panel shows the same after including 
feedback. 
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Fig. 21.— 



Sketch of the accretion disk-corona system as derived in this paper. 
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TABLE 1 
Parameters of the Light Curves. 





Data set 


Start 


End 


T (days) 1 


AT (days) 2 


^points 


X-ray 


Longlook 
Medium 
Monitor 


2002 December 13 
2006 November 
2002 March 


2002 December 14 
2007 January 
2007 May 


1.5 
60.0 
1910.0 


0.01 
0.25 
15.0 


150 
240 
1050 


Optical 


Monitor 


2004 August 


2008 January 


1250.0 




154 


Radio 


Monitor 


2002 March 


2008 January 


2167.0 
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"Total length of light curves 
'Bin size 



TABLE 2 

Time, Area and Width of the X-ray Dips and 37 GHz Flares, and Times of Superluminal Ejections. 



Parameters of X- 
Timc (start) Time (min.) 


ray Dips 
Area 1 


Width 2 


Ejection Times 
To 


Knot ID 


Parameters 
Time (peak) 


of 37 GHz 

Area 3 


Flares 
Width 2 


_4 


2002.15 


_4 


_4 


2002.23 ± 0.03 


02A 


2002.39 


195 


95 


2002.19 


2002.30 


657 


120. 


2002.65 ± 0.04 


02B 


2002.58 


404 


75 


2002.75 


2002.76 


12 


5. 












2003.02 


2003.12 


29 


22.5 


2003.35 ± 0.15 


03A 


2003.35 


17 


17 


2003.32 


2003.58 


172 


52.5 


2003.67 ± 0.02 


03B 


2003.72 


441 


85 


2003.66 


2003.82 


184 


40. 


2003.81 ± 0.03 


03C 


2003.92 


24 


20 


2003.95 


2003.98 


44 


15. 


2003.98 ± 0.03 


03D 








2004.12 


2004.17 


20 


15. 


2004.16 ± 0.05 


04A 








2004.21 


2004.37 


364 


65. 


2004.37 ± 0.03 


04B 


2004.38 


126 


52 














2004.49 


60 


25 


2004.62 


2004.66 


30 


12.5 


2004.82 ± 0.05 


04C 


2005.05 


61 


27 


2005.09 


2005.12 


33 


17.5 


2005.14 ± 0.03 


05A 


2005.23 


21 


15 


2005.19 


2005.39 


351 


65. 


2005.34 ± 0.02 


05B 


2005.36 


10 


10 














2005.57 


92 


57 














2005.80 


127 


57 


2005.94 


2006.04 


287 


70. 


2006.00 ± 0.03 


06A 


2006.43 


126 


105 


2006.38 


2006.44 


382 


72.5 


2006.72 ± 0.05 


06B 


2006.88 


598 


62 


2006.96 


2007.09 


193 


62.5 


2007.05 ± 0.02 


07A 


2007.43 


132 


72 










_4 




2007.87 


122 


60 



"Equivalent width of dips. Units: 10 — 6 erg cm~ 2 
'Units: days 

c Area under the flare light curves. Units: 10 — 6 erg cm -2 
^Insufficient data 



